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Abstract 

The  aim  of  this  study  is  the  understanding  of  the  main  phenomena  governing  fuel  cell  performances.  We  present  a  fuel  cell  model  that 
takes  into  account  gas  diffusion  in  the  porous  electrodes,  water  diffusion  and  electro-osmotic  transport  through  the  polymeric  membrane,  and 
heat  transfer  in  both  the  Membrane  Electrodes  Assembly  (MEA)  and  bipolar  plates.  This  model  is  constructed  by  combining  independent 
descriptions  of  heat  and  mass  transfers  in  the  cell  with  a  third  description  of  coupled  charge  and  mass  transfers  in  the  electrodes,  considered 
porous. 

The  results  show  that  thermal  gradients  in  the  MEA  could  lead  to  thermal  stresses  at  high  current  densities.  The  feeding  gas  temperature 
influence  on  the  cell  temperature  is  also  important. 

©  2005  Elsevier  B.Y.  All  rights  reserved. 
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1.  Introduction 

Compared  to  thermodynamic  machines  used  for  electric¬ 
ity  production,  one  of  the  main  advantages  associated  with 
fuel  cells  is  that  they  are  not  limited  by  Carnot  efficiency. 
Besides,  no  moving  part  is  required  to  convert  chemical 
energy  into  electric  energy.  In  this  paper,  we  consider  Proton 
Exchange  Membrane  Fuel  Cells  (PEMFCs):  PEMFCs  are 
composed  of  a  polymer  electrolyte  sandwiched  between 
two  porous  backing  layers  to  form  a  Membrane  Electrode 
Assembly  (MEA).  The  electrodes  (the  active  layers)  are 
inserted  between  the  electrolyte  and  the  backing  layers. 
Following  Springer  et  al.  [1],  we  refer  to  the  backing  and 
active  layers  (considered  together  but  without  the  membrane) 
as  the  Gas  Diffusion  Electrodes  (GDE).  The  MEA  is  located 
between  two  graphite  bipolar  plates,  which  collect  electric 
current  and  allow  gas  feeding  and  system  cooling.  At  the 
anode,  EL  is  oxidized  (1),  liberating  electrons  and  producing 
protons.  The  electrons  flow  to  the  cathode  via  an  external 
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circuit,  where  they  combine  with  the  protons  and  the  oxidant 
(oxygen)  to  produce  water  (2).  Proton  transfer  from  the 
anode  to  the  cathode  through  the  polymer  membrane  closes 
the  electrical  circuit.  The  overall  reaction  (3)  differs  from 
hydrogen  combustion  because  part  of  the  Gibbs  free  energy 
can  be  converted  into  electric  energy 


anode 

H2  —  2H+  +  2e~ 

(1) 

cathode 

\<d2  +  2H+  +  2e_  - 

-*  h2o 

(2) 

overall 

reaction  H2  +  Ao2  - 

*  h2o 

(3) 

A  careful  water  management  is  necessary  to  ensure  that 
the  membrane  remains  fully  hydrated,  in  order  to  improve 
the  ionic  conductivity  and  to  avoid  electrodes  flooding.  This 
occurs  when  an  excess  of  liquid  water  restrains  the  species 
access  to  the  active  layers.  Water  saturation  pressure  being 
a  strong  and  non-linear  function  of  temperature,  water  man¬ 
agement  is  closely  linked  to  thermal  management.  Both  water 
and  thermal  managements  are  key  issues  for  PEMFCs  sys¬ 
tems. 

Various  models  help  to  understand  the  MEA  behaviour. 
According  to  Bernardi  and  Verbrugge  [2]  and  Jaouen  et  al. 
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Nomenclature 

bi 

tafel  slope  (V  dec-1) 

cn2o 

water  concentration  (molm3) 

^pore 

pore  diameter  (m) 

Dij 

diffusion  coefficients  of  species  i  to  j  (m2  s-1) 

Dm 

effective  diffusion  coefficients  of  water  in  the 
membrane  (m2  s-1) 

Du 

effective  diffusion  coefficients  (m2  s-1) 

Dfj 

Knudsen  diffusion  coefficient  (m2  s-1) 

EW 

equivalent  weight  (kg  mol-1) 

T^cell 

cell  potential  (V) 

E° 

standard  potential  (V) 

F 

Faraday’s  constant  (C  mol-1) 

h 

convective  heat  transfer  coefficient 
(W  m-2  K-1) 

HRa,c 

relative  humidities  (-) 

if 

Faradic  current  density  (A  m-2) 

*o(0 

exchange  current  density  (A  m-2) 

Tcell 

total  current  density  (A  m-2) 

L 

thickness  (m) 

L^o 

latent  heat  of  water  evaporation  (kJ  mol-1) 

Mi i 

molar  weight  (kg  mol-1) 

N 

density  of  molar  flux  (mol  m-2  s) 

Nu 

Nombre  de  Nusselt  (-) 

P 

pressure  (Pa) 

Pi 

partial  pressure  (Pa) 

P  sat 

saturated  vapor  pressure  (Pa) 

Q 

surfacic  heat  sources  (W  m-2) 

Qj 

volumetric  heat  source  in  the  membrane 
(Wm"3) 

R 

Universal  gas  constant  (J  mol-1  K-1) 

yi 

molar  rate  of  the  gas  species  i  (-) 

z 

thermal  resistance  (W  m-2) 

Greek  letters 

y 

roughness  factor  (-) 

8 

size  of  the  agglomerate  (m) 

£ 

porosity  (-) 

ceff 

&a,c 

effective  porosity  (in  the  agglomerate)  (-) 

m 

overpotential  (V) 

K 

thermal  conductivities  (W  m-1  K-1) 

X 

water  content  in  the  membrane  (-) 

1 

tortuosity  (-) 

Pdry 

dry  Nation®  density  (kgm-3) 

°H+ 

protonic  conductivity  (S  m-1) 

<Teff 

au+ 

effective  protonic  conductivity  (S  m-1) 

Jo 

osmotic  coefficient  (mol  mol-1) 

0 

thermal  flux  (W  m  2) 

Subscripts 

a 

anode 

act 

active  layers 

b 

backing  layers 

BP 

bipolar  plates 

c 

cathode 

h2 

hydrogen 

h2o 

water 

m 

membrane 

n2 

nitrogen 

o2 

oxygen 

[3],  since  water  is  produced  at  the  cathode,  the  impact  of 
water  management  on  cathode  performances  is  significant  at 
high  current  densities.  However,  these  authors  do  not  model 
mass  transfer  in  the  whole  MEA.  Experimental  works  of 
Buchi  and  Srinivasan  [4]  confirmed  the  importance  of  gas 
hydration  and  Sena  et  al.  [5]  emphasized  the  links  between 
water  transport  in  the  membrane,  membrane  thickness  and 
feeding  gases  composition.  According  to  these  authors,  the 
ionic  transport  through  the  membrane  is  the  limiting  phe¬ 
nomenon  when  the  gas  hydration  is  low.  On  the  other  hand, 
they  also  showed  that  the  access  of  gases  to  catalytic  sur¬ 
faces  becomes  a  limiting  factor  when  the  membrane  is  suffi¬ 
ciently  hydrated  [5].  Costamagna  [6,7]  and  Okada  et  al.  [8], 
who  modelled  the  charge  and  mass  transfers  inside  the  MEA 
in  steady  state  or  in  transient  conditions,  confirmed  these 
observations. 

The  literature  proposes  various  ways  to  predict  the 
electrochemical  behaviour  of  electrodes.  Springer  et  al. 
[1,9]  proposed  analytical  formulations  and  emphasized  the 
capacitive  and  resistive  effects  of  oxygen  reduction.  Song 
et  al.  [10]  and  Eikerling  and  Kornyshev  [11]  completed 
these  information  by  their  descriptions  of  the  electrode 
microstructure.  This  must  be  considered  in  order  to  explain 
the  electrodes  electrochemical  behaviour  and  to  estimate 
the  overvoltage.  Thus,  Rowe  and  Li  [12]  took  into  account 
the  volumetric  characteristics  of  the  active  layer  while 
considering  it  as  a  homogeneous  and  non-porous  mixture  of 
polymer,  carbon  and  catalyst  [13]. 

The  model  of  Rowe  and  Li  [12]  also  predicts  the  tempera¬ 
ture  gradients  inside  the  MEA  as  well  as  their  consequences 
on  PEMFC  performances.  Similarly,  Fuller  and  Newman 
[14]  showed  that  temperature  distributions  depend  on  mass 
balance.  They  emphasized  the  need  to  control  well  the 
membrane  hydration  and  PEMFC  operating  temperature. 
Yan  et  al.  [15]  show  that  the  temperature  could  affect  water 
distribution  inside  the  membrane.  This  can  be  damaged  at 
high  current  density. 

Although  most  of  the  models  used  are  one-dimensional, 
they  correctly  predict  the  electrochemical  behaviour  of  MEA. 
However,  Costamagna  [6,7]  and  Djilali  and  Lu  [16]  suggested 
that  a  multidimensional  model  could  improve  the  description, 
since  gas  composition  and  temperature  vary  along  the  feeding 
channels. 
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The  model  presented  in  this  work  describes,  in  steady 
state,  one-dimensional  heat  and  mass  transfers  in  the  whole 
cell  and  charge  and  mass  transfers  in  the  electrodes.  It  is  con¬ 
structed  by  combining  independent  descriptions  of  heat  and 
mass  transfers  in  the  cell  with  a  third  description  of  coupled 
charge  and  mass  transfers  in  the  electrodes.  Results  of  the 
whole  cell  mass  transfer  model  are  used  as  boundary  con¬ 
ditions  in  the  electrode  model,  which  provides  information 
about  heat  sources,  necessary  to  determine  heat  fluxes  in  the 
MEA.  The  heat  transfer  model  also  requires  data  about  wa¬ 
ter  fluxes  or  gas  velocity  coming  from  the  cell  mass  trans¬ 
fer  model.  Finally,  this  description  of  heat,  mass  and  charge 
transfers  allows  the  prediction  of  the  thermal,  concentration 
and  potential  fields  in  a  single  cell. 

2.  Theoretical  models 

2.7.  Geometric  description 


MEA,  the  electronic  and  ionic  transfers  occurring  in  these 
areas  are  assumed  to  be  without  significant  consequence  on 
the  mass  transfer  in  the  whole  cell.  Charge  and  mass  trans¬ 
fers  in  the  electrodes  are  considered  in  order  to  evaluate  their 
electrochemical  behaviour.  The  overvoltages  are  then  inter¬ 
preted  in  term  of  heat  sources  in  the  heat  transfer  model.  Let 
us  note  that  mass  and  charge  transfers  are  always  calculated 
assuming  that  the  MEA  is  isothermal.  The  operating  condi¬ 
tions  are  always  chosen  in  order  to  avoid  water  condensation 
in  the  backing  layers  and  the  gas  mixtures  are  assumed  to 
obey  the  ideal  gas  law. 

2.2.  Mass  transfer  in  the  whole  MEA 

The  source  terms  are  the  molar  flux  densities  of  hydrogen 
and  oxygen  consumed  (negative  sign)  at  the  anode  and  at 
the  cathode,  respectively,  and  the  molar  flux  density  of  water 
produced  at  the  cathode  (positive  sign).  They  are  given  by 
the  following  expressions: 


A  PEMFC  single  cell  can  be  described  schematically 

(Fig.  1)  as  an  assembly  of  seven  layers  constituting  four  dis¬ 
tinct  areas: 

•  On  both  sides  of  the  MEA,  the  bipolar  plates  cumulate 
three  functions:  thermal  control  thanks  to  cooling  water 
channels,  gas  feeding  via  the  gas  channels  and  current  col¬ 
lector. 

•  The  porous  backing  layers  made  of  carbon  fibres  coated 
with  hydrophobic  PTFE  ensure  a  homogeneous  gas  distri¬ 
bution  on  the  electrodes  and  limit  liquid  water  accumula¬ 
tion  in  the  MEA. 

•  The  active  layers  are  thin  layers  of  platinum  dispersed  on 
carbon  particles  embedded  in  the  polymer  electrolyte.  The 
catalyzed  electrochemical  and  chemical  reactions  occur  in 
these  layers. 

•  Only  protons  and  water  can  go  through  the  polymer  mem¬ 
brane  located  at  the  centre  of  the  assembly. 


The  membrane  being  impervious  to  gases,  hydrogen  and 
oxygen  only  flows  from  the  bipolar  plates  to  the  anode  and 
to  the  cathode,  respectively.  Only  water  can  go  through  the 
membrane,  in  a  direction  depending  on  current  density  and 
on  gas  hydration  (5)  [8].  Fluxes  are  considered  positive  in  the 
x-direction  (Fig.  1) 

^h2o  =  +  ^h2o  and  ^h2o  =  ^h2o  (5) 

Most  of  the  models  encountered  in  the  literature  use 
Stefan-Maxwell  equations  to  describe  the  diffusion  of  gas 
in  the  backing  layers  (6)  [1]: 


dy; 

dx 


yiNj  -  yjNf 


Considering  the  low  thickness  of  the  active  layers,  more 
than  100  times  thinner  than  the  others  components  of  the 


Membrane-Electrode  Assembly 

Anodic  Gas  Cathodic  Gas 

Diffusion  Electrode  Diffusion  Electrode 


/" -  - >  ^ -  - "N 

Anodic  Cathodic  02 


1 - 1 - 1 - !—►  x 

0;i  xc  Lc  Om  xm  Lm  Oc  xc  Lc 


Fig.  1.  Schematic  description  of  the  PEMFC  single  cell. 


and 


yi  =  1  at  both  anode  and  cathode  sides  (7) 

i 


with  /,  j  —  H2,  H^Ogas  at  the  anode  side  and  /,  j  —  O2,  N2, 
H^Ogas  at  the  cathode  side. 

The  backing  layers  being  porous,  it  is  necessary  to  use 
effective  diffusion  coefficients  (8)  [17]: 


Defl  = 

IJ 


2/3 


The  anode  and  cathode  thickness  are  neglected:  the  elec¬ 
trodes  are  considered  as  interfaces  between  the  membrane 
and  the  backing  layers. 

Considering  the  molar  flux  densities  of  water  in  the  back¬ 
ing  layer  (4)  as  parameters,  and  using  the  concentrations  of 
gases  in  the  bipolar  plates  feeding  channels  as  boundary  con¬ 
ditions,  Eq.  (6)  lead  to  the  following  analytic  expressions  for 
the  species  concentration  distributions  (9)-(14). 
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In  the  anodic  backing  layer: 
if  jVh2  +  A^H2Oa  7^  0 


;m2(xc)  =  yN2 


exp(fcnxc) 

exp(knl) 


(13) 


jH2fe)  = 


Afir 


Afit2  +  Nu2o, 


+  Ah, 


Afir 


Nh2  +  A^h2o£ 


exp(kaxa) 

(9) 


with 


[A^H2]  +  A^H2oa 


cD 


eff 

h2,h2o 


5 


p 

c  =  — ,  }’H2Oa(A'a)  =  1  -  )’H2  lA'a) 


if  jVh2  +  A'W).,  =  0 


jHjCZa)  = 


Nhz 


0 


rDeff 

c^h2,h2o 


x&  + 


with  yH2OaUa)  =  1  -  V'H2(-ta). 

In  the  cathodic  backing  layer: 


if  kn  -  kc  /  0 


yN2(*c)  =  )>n2 


exp(fcnxc) 

exp(£n/) 


with 


Afit2oc 


cD 


eff 

n2,h2o 


+ 


cD 


eff 

n2,o2 


(10) 


(11) 


yo2(xc )  = 


cD 


eff 

o2,h2o 


y^2NQi  (  1 _ 1  exp(kcxc) 

kn  ~  kc  yC^o2,H20  cDoIn2  J  J  exP(kcO 


+ 


1 


cD 


eff 

o2,h2o 


1 


rDeff 

cuo2^2 


exp(Anvc) 

exp(An/) 


+ 


cD 


eff 

o2,h2o 


1 

kc 


(12) 


with 


Nh2Oc  +  Nq2 


cD 


eff 

o2,h2o 


and 


yu2 oUc)  =  1  -  X>2C*c)  -  W2(^c) 


if  kn  —  kc  =  0 


w2(*c)  =\y°o2-  j-.eff 


Nq2  1  \  exp(£cvc) 


cDo2,h2o  kc  J  exP(^cO 


.0 


+  y^oN0: 


1 


1 


rDeff 

c/Ad2,h2o 


rDeff 

CLfo2,N2 


X 


xc  exp(Acvc)  -  /  exp(Ac/) 
exp(£n/) 


+ 


NO: 


r  Deff  k 

cVo2,u2o  Kc 


1 


(14) 


where  yj^ ,  }?q2  and  are  molar  concentrations  in  the  gas 
channels. 

Bernardi  and  Verbruge  [2]  describe  water  transport  in  the 
polymer  membrane  using  Darcy  equations  completed  by  a 
supplementary  term  representing  the  electro-osmotic  water 
transfer  due  to  electric  field  and  proton  flux.  We  do  not  retain 
this  description,  which  needs  a  high  number  of  parameters, 
some  of  them  difficult  to  evaluate.  The  phenomenological 
model  of  water  transport  in  the  membrane  proposed  by 
Springer  et  al.  [1]  is  used  here.  In  this  model,  the  water 
content  (A)  of  the  membrane  is  defined  as  the  number  of 
water  moles  per  mole  of  sulfonic  acid 

EW 

A  = - cr2o  (15) 

Pdry 

where  EW  (equivalent  weight)  represents  the  dry  membrane 
weight  per  mole  of  sulfonate  group  (kg  mol-1);  pdry  is  the 
density  of  dry  polymer  (kg  m~3). 

According  to  Okada  et  al.  [8],  the  water  electro-osmotic 
flux  through  the  membrane,  always  directed  from  anode  to 
cathode,  is  a  linear  function  of  the  proton  flux  imposed  by 
the  current  density: 

^osmotic  =  T  ^11  and  T  =  lro  (16) 

where  r  stands  for  the  water  molecules  number  dragged  by  a 
proton  [8].  rG  is  the  electro-osmotic  drag  coefficient  (Table  1). 


Table  1 

Parameters  used  in  mass  transfer  model 


Parameter 

Value 

Refs. 

A  a  =  Lc 

230  x  10-6  m 

This  work 

Lm 

125  x  10“6m 

This  work 

£back 

0.8 

[1] 

EW 

1 . 1  kg  mol- 1 

[8] 

Pdry 

2020  kg  m-3 

[8] 

To 

2.5/22 

[8] 

Dm 

3  x  10“ 10  m2  s-1 

[8] 

Deff 

^h2/h2o 

1.63  x  10-4  m2  s-1 

[20] 

Deff 

^o2/h2o 

3.20  x  10-5  m2  s-1 

[20] 

Deff 

^o2/n2 

2.41  x  10-5  m2  s-1 

[20] 

Deff 

^h2o  /  n2 

3.35  x  10“5  m2  s_1 

[20] 

420 


J.  Ramousse  et  al.  /  Journal  of  Power  Sources  145  (2005)  416-427 


The  water  diffusion  flux  through  the  membrane  de¬ 
pends  linearly  on  the  water  concentration  gradient.  The  sum 
of  diffusion  and  electro-osmotic  fluxes  yields  a  first  order 
differential  equation  in  the  A  variable: 


Nm 

7Vh2o 


a  r  osmotic  i  Ardiff 
7VH20  "h  7VH20 


Tcell 


^  Pdry  dA 

ly  rn 

EW  dx 


(17) 


And  the  solution  of  Eq.  (17)  is: 


Nm  — 

7Vh2o  - 


Jo  J cell 

F 


Aa  + 


A  r  —  A 


a 


1  exp(AmLm) 


(18) 


with 


EW  TqJqqW 


Pdry  Dm  F 


Ac  and  Aa  stand  for  polymer  water  content  at  the  mem¬ 
brane/cathode  and  membrane/anode  interfaces,  respectively. 
Then,  assuming  thermodynamic  equilibrium  between  water 
vapor  in  the  backing  layers  and  liquid  water  in  the  polymer, 
the  sorption  curve  of  Hinatsu  et  al.  [18]  is  used  to  associate 
partial  pressure  values  to  Aa  and  Ac: 


A  =  0.3  +  10.8 


2 

+  14.1 


3 


(19) 


The  solving  algorithm  determines  the  water  concentra¬ 
tions  at  the  GDE/membrane  interfaces  that  satisfy  Eq.  (5). 
The  membrane  ionic  conductivity  depends  on  its  temperature 
and  water  content,  according  to  the  correlation  of  Neubrand 

[19]: 

<th+  =  e("£A(1/;_1/353))(0.00  13A3  +  0.0298A2  +  0.2658A) 

(20) 


possible  to  investigate  in  details  the  physicochemical  mech¬ 
anisms  taking  place  in  the  electrodes: 


•  Gas  concentration  variation  in  the  pores. 

•  Agglomerate  protonic  conductivity  dependence  on  hydra¬ 
tion. 

•  Microstructural  effects  (platinum  loading,  agglomerate 
size). 

However,  there  are  phenomena  disparities  linked  to  water 
transfers  in  the  fuel  cell,  since  under  high  current  densities  the 
anode  will  tend  to  be  dehydrated  and  cathode  to  be  flooded. 
Under  these  conditions  the  model  assumptions  are  not  iden¬ 
tical  for  both  electrodes. 

Fig.  2  represents  schematically  the  simplified  pore  struc¬ 
ture  assumed  in  the  agglomerate  model.  The  solid  phase 
consists  of  a  homogeneous  mixture  of  polymer  electrolyte 
(Nation®),  carbon,  and  catalyst  (Pt). 

The  active  layer  thickness  Lact,  the  agglomerate  (assumed 
cylindrical)  diameter  8  and  the  porosity  £ac t  characterize  the 
active  layer  geometry.  The  active  layer  porosity  £act  is  defined 
by: 


£act  — 


pore 


’pore 


Epore  T  E 


agg 


’pore 


+  5 


agg 


(23) 


where  Sagg  is  the  agglomerate  section  and  .Spore  is  the  pore 
section. 

And  the  pore  equivalent  diameter  can  be  estimated  by: 


dp  ore  =  -2^—S  (24) 

As  mentioned  by  Middelman  [24],  the  best  electrodes 
performances  are  achieved  with  the  highest  possible  ratio 
Lact/3.  The  electrodes  electrochemical  behaviour  also 


with 

Ea  =  2640  e(_a6A)  +  1183  (21) 


And  the  ohmic  drop  is  given  by: 

fLm  1 

0m  =  Tcell  /  ;  7  dx  (22) 

JO  °H+  (*) 

The  values  of  the  main  parameters  used  in  this  model  are 
summarized  in  Table  1 . 

2.3.  Active  layers  equations 

Considering  the  results  of  a  cylindrical  model  of  pores  in¬ 
side  the  electrode,  Srinivasan  and  Hurwirtz  [21]  assert  that 
it  is  preferable  to  use  a  volumetric  description  of  the  active 
layer  and  a  formalism  of  the  Butler- Volmer  type  than  to  use 
a  description  derived  from  Tafel  law.  The  agglomerate  model 
[22]  that  assumes  the  presence  of  macro-pores  in  the  elec¬ 
trodes,  seems  to  be  well  adapted  to  the  PEMFC  anode  and 
cathode.  Siegel  et  al.  [23]  used  it  to  simulate  the  electrochem¬ 
ical  behaviour  of  a  MEA.  Thanks  to  this  kind  of  model,  it  is 
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Fig.  2.  Schematic  description  of  the  active  layers. 
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depends  on  a  last  parameter,  the  specific  area,  which  is  a 
function  of  the  catalyst  load  Wpt  (g  m-2)  and  is  recurrently 
expressed  as  follows: 


active  platinum  area 
MEA  section 


(25) 


In  the  agglomerate,  both  an  ionic  current  density  and 
an  electronic  current  density  flow.  They  vary  with  the  v 
coordinate.  The  ionic  current  density  is  nil  at  the  active 
layer/backing  layer  interfaces  (x  =  0)  and  the  electronic  den¬ 
sity  is  nil  at  the  membrane/active  layers  interfaces  ( x  =  Lact ). 
In  steady  state,  the  local  variation  of  the  electronic  current 
density  is  the  opposite  of  the  local  variation  of  the  ionic  cur¬ 
rent  density;  it  is  also  equal  to  the  faradaic  current  density  if. 
The  cell  current  density  /ceii  is  the  sum  (over  the  electrode 
thickness  Lact)  of  the  faradaic  current  density: 


(26) 


2.3.1.  Anode  model 

PEMFCs  anode  is  seldom  studied.  Boyer  et  al.  [25]  assert 
that  the  ionic  conductivity  of  the  active  layer  influences  the 
electrochemical  performances  of  the  electrodes.  Further¬ 
more,  due  to  electro-osmotic  water  flux  in  the  membrane, 
the  polymer  of  the  anode  agglomerate  is  more  prone  to 
dehydration  than  that  of  the  cathode  [26].  Considering  this, 
it  seems  appropriate  that  the  anode  model  takes  into  account 
both  the  ohmic  drop  in  the  active  layer  and  the  hydrogen 
concentration  variations  in  the  gas  pores.  This  modelling  of 
the  ionic  transfer  completes  the  description  of  Rowe  and  Li 
[12]  for  mass  transfer  in  the  non-porous  active  layer. 

The  main  assumptions  of  the  anode  model  are  the  follow¬ 
ing: 

1 .  Mass  transfer  of  H2O  and  H2  is  considered  only  in  the  gas 
pores,  there  is  no  diffusion  of  water  in  the  agglomerate. 

2.  The  H2O  concentration  in  the  agglomerate  is  imposed  by 
thermodynamic  equilibrium  with  the  gas  phase. 

3.  The  ionic  conductivity  of  the  agglomerate  is  a  function  of 
the  ionic  conductivity  of  the  polymer,  which  depends  on 
its  water  content. 

4.  Only  the  migration  of  protons  is  considered  in  the  agglom¬ 
erate. 

5.  Pressure  and  temperature  are  uniform. 

The  Stefan-Maxwell  equations  [1]  completed  by  a  sup¬ 
plementary  term  for  Knudsen  diffusion  in  the  small  pores  of 
the  active  layer  [17]  describe  diffusion  of  hydrogen  and  water 
in  the  gas  phase: 

d yi  _  RT  ( y[N j  —  yjNi 
d x~~P~  \^j  [  ~Dij 


the  Knudsen  diffusion  coefficient  D\  being  expressed  as  fol¬ 
lows  [17]: 


Dk  —  d 

LJi  —  ^pore 


I 

5^act 


(28) 


with  dpore  and  £act  given  by  Eqs.  (23)  and  (24). 

The  hydrogen  consumption  in  the  gas  pores  of  the  active 
layer  yields  the  following  mass  balance  equation: 


9/vh2  y  if 
dx  Lact  2 F 


And  the  water  mass  balance  leads  to: 


(29) 


3*h2o 

dx 


=  0 


(30) 


The  boundary  conditions  of  Eqs.  (27),  (29)  and  (30)  are 
given  by  the  output  of  the  mass  transfer  model  in  the  whole 
cell: 


-  Concentration  of  each  species  at  the  active  layer/backing 
layer  interfaces  (9)-(10). 

-  Molar  flux  density  of  water  and  hydrogen  (4)-(5). 


The  overpotential  distribution  Eq.  (33)  in  the  active  layer 
can  be  calculated  by  means  of  the  proton  balance  (31)  and 
Ohm’s  law  (32)  in  the  ionic  phase: 


di 

dx 


r 


L 


H 


act 


dr] a  _  J_ 
dx  ~  crj'ff 


(31) 

(32) 


(33) 


The  effective  ionic  conductivity  of  the  agglomerate  de¬ 
pends  on  Nation® ’s  conductivity  [19]  (20)  and  (21),  which 
is  simply  corrected  as  a  function  of  the  agglomerate  polymer 
content,  £Nafion: 


^H+W  =  SNafionOH+M 


(34) 


Since  the  polymer  water  content  A  depends,  via  the  sorp¬ 
tion  curves  (19),  on  water  partial  pressure  in  the  gas  pores, 
the  agglomerate  ionic  conductivity  is  thus  a  function  of  the  v 
co-ordinate;  differential  equation  (31)  has  to  be  solved  with 
the  following  boundary  conditions: 


The  ionic  current  density  is  nil  at  the  active  layer/backing 


layer  interface  cr^+  ^ 


=  0 


Electrode  overpotential  at  the  membrane/active  layer  in¬ 
terface  is  known  ( r](x  =  Lact)  =  ijo)-  In  practice,  rjo  is  ad¬ 
justed  until  /Ceii  equals  the  value  used  for  the  determination 
of  hydrogen  consumption  (4). 


The  solution  of  Eqs.  (31)— (33)  yields  the  distributions  of 
concentrations  and  electrodes  overpotential  as  function  of  x- 
coordinate. 


422 


J.  Ramousse  et  al.  /  Journal  of  Power  Sources  145  (2005)  416-427 


Two  mechanisms  are  commonly  associated  with  hydrogen 
oxidation  on  platinum  in  acid  mixture:  Heyrovsky-Volmer 
and  Tafel-Volmer  [27].  With  a  fast  Volmer  step,  which  is 
usually  the  case  in  strong  acid  media,  the  oxidation  kinetics 
can  be  written  as  follows: 

if  =  /0(H2)  (e(2-3^W^  _  e-(2.3 /**)*« 

V  yH2 

where  yj^?  is  the  hydrogen  concentration  at  the  open-circuit. 
2.3.2.  Cathode  model 

According  to  the  majority  of  the  authors,  fuel  cells  lack 
of  efficiency  (in  terms  of  electrochemical  conversion)  is  due 
mostly  to  phenomena  occurring  at  the  cathode.  Kim  et  al. 
[28]  compared  the  results  of  their  models  to  experimental 
data  in  various  operating  conditions  and  concluded  that  the 
oxygen  reduction  process  is  limiting.  The  cathode  being 
the  place  of  water  production,  it  is  reasonable  to  assume 
that  ionic  resistivity  of  the  well-hydrated  polymer  can  be 
neglected.  However,  oxygen  access  to  the  catalyst  sites 
remains  problematic,  not  only  because  of  the  flooding 
risk,  but  also  due  to  its  low  partial  pressure.  As  Bernardi 
and  Verbruge  [2]  proposed,  a  one-dimensional  model  can 
describe  the  cathode  operation  by  taking  into  account  only 
oxygen  diffusion  in  the  gas  pores  and  reduction  kinetics. 

The  main  assumptions  of  the  cathode  model  are  the  fol¬ 
lowing: 

1 .  Mass  transfer  of  H2O  and  O2  is  considered  only  in  the  gas 
pores. 

2.  Diffusion  of  water  in  the  agglomerate  is  neglected. 

3.  The  ohmic  drop  in  the  active  layer  is  neglected. 

4.  Pressure  and  temperature  are  uniform. 

Gas  diffusion  inside  the  cathodic  gas  pores  is  described 
by  Stefan-Maxwell  equations  modified  to  take  account  of 
Knudsen  diffusion  ((27)  and  (28)).  Oxygen  consumption  and 
water  production  in  the  porous  active  layer  yield  the  following 
mass  balance  equations: 


Table  2 


Parameter  used  in  the  active  layer  models 


Parameter 

Value 

Refs. 

L 

10  x  10-6  m 

This  work 

<5 

3  x  10-6  m 

[23] 

£act 

0.3  m3  m~3 

This  work 

^Nafion 

0.5  m3  m-3 

This  work 

1.7 

This  work 

y 

100  m2  m“2 

[27] 

bn2  — 

30  mV  dec-1 

[27] 

*0(H2) 

1  x  10~2  A  cm-2 

[27] 

bo2  =  b*0i 

120  mV  dec-1 

[27] 

io(02) 

4  x  10-7  A  cm-2 

[27] 

A  Butler- Volmer  law  [3,22,29]  is  generally  used  to  de¬ 
scribe  oxygen  reduction  kinetics  in  acid  media  (Nation®): 


k(2.3/6o9)|?7c|  yOi  _  -(2.3/fc*  Jli/c 


if  =  fl)(02)  |  e— ' -  e  — o2-/ci 

yo: 


(39) 


2.3.3.  Main  parameters  and  overvoltage  calculation 

The  cell  potential  is  determined  by  means  of  the  following 
expressions: 

^cell  —  E  —  ?7a  U=0  "T  7c  —  7  m  (40) 


A G(T,  P ) 
IF 


(41) 


where  A G(T,  P )  is  given  by  the  Nernst  equation  and  rjm  is 
a  result  of  the  cell  mass  transfer  model  (22).  r{ ja  (v  =  0)  (35) 
and  T)c  (39)  are  estimated  thanks  to  the  GDE  models  with 
parameters  whose  values  are  reported  in  Table  2. 


2.4.  Heat  transfer  in  the  single  cell 

Various  phenomena  are  responsible  for  heat  production  by 
fuel  cells  [30] .  Firstly,  the  ionic  resistivity  of  the  membrane 
through  which  protons  flow  is  at  the  origin  of  a  volumetric 
heat  source  distributed  in  its  whole  thickness  (22): 


Qj  —  *^cell7m 


(42) 


dx 


M02  |  y  ^ 

dx  Lact  4  F 


Mh2o 

dx 


y  H 


Fact  ^ F 


=  0 


(36) 

(37) 

(38) 


Once  again,  the  boundary  conditions  of  Eqs.  (27), 
(36)-(38)  are  given  by  results  of  the  mass  transfer  model 
in  the  whole  cell: 


Secondly,  electrochemical  reactions,  and  water  sorption  or 
desorption,  imply  heat  production  (or  even  consumption  in 
some  particular  cases)  at  the  electrodes.  The  heat  production 
mechanisms  occurring  at  the  electrodes  are  the  following: 

-  Global  reaction  entropy.  Even  in  fully  reversible  condi¬ 
tions,  the  simultaneous  hydrogen  oxidation  at  the  anode 
and  oxygen  reduction  at  the  cathode  produce  heat  (43) 

h2  +  ±o2  -*  H2Ov, 

TAS  =  —4.7  kJ  mol-1  at  T  =  298  K 


-  Concentration  of  each  species  at  the  backing  layer/active 
layer  interfaces  (1 1)— (14). 

-  Molar  flux  density  of  water  and  oxygen  (4)-(5). 


H2  +  ±02  H2Oi, 

TAS  =  -48.7  kJ  mol-1  at  T  =  298  K 


(43) 
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/Sreac  ,  /Sreac  _  ^cell  TA  S 

*2a  '  ^ c  2^7 

It  is  not  easy  to  determine  which  part  of  the  global  reaction 
entropy  must  be  attributed  to  each  electrode:  a  quick  review 
of  the  literature  shows  that,  according  to  the  authors,  the 
half-reaction  of  hydrogen  oxidation  (1)  can  be  isothermal 
(0aeac  =  0)  [31],  endothermic  (0aeac  <  0)  [32]  or  exother¬ 
mic  (0fac  >  0)  [33].  In  this  paper  the  hydrogen  oxidation 
is  assumed  isothermal  and,  consequently,  the  entropy  of 
oxygen  reduction  (2)  equals  that  of  water  formation  reac¬ 
tion  (43). 

-  Electrode  overpotentials  result  from  the  electrodes  models 
((35)  and  (39)).  The  associated  heat  production  is  given  by 
Eq.  (44).  The  thickness  of  the  electrodes  is  neglected  in  the 
heat  transfer  model 

G°VerP0‘  =  >h(*  =  0)  ■/ cell 


g  overpot  =  lVcUcell  (44) 

-  Water  sorption  or  desorption  at  the  GDE/membrane  inter¬ 
faces 


e|orp  =  A//sH2°Af2o  and 

gsorp  =  _  A Hf2°N^l0  with  A  h's'-°  >  0  (45) 


-  Water  being  liquid  in  the  polymer  membrane  and  gaseous  in 
the  backing  layers,  the  sorption  and  desorption  phenomena 
are  responsible  for  heat  production  or  absorption,  respec¬ 
tively.  The  enthalpy  of  sorption  A  H^2°  is  assumed  equal 
to  the  latent  heat  of  water  vaporization  ALy2°  [34]. 


In  steady  state,  heat  transfer  through  the  MEA  is  modelled 
using  the  energy  equation  applied  to  incompressible  fluids 
with  constant  physical  properties  in  porous  media  [35].  The 
volumetric  heat  source  Qj  is  nil  in  the  backing  layers 


d2T 


K 


+  ( Qj ) 


(46) 


where  Si  is  the  i  species  volumetric  fraction  in  the  backing 
layer,  defined  such  as  ^  =  £b,m-  k  is  the  equivalent  ther¬ 

mal  conductivity  of  the  porous  media  (Table  3).  The  molar 


Table  3 

Parameter  used  in  heat  transfer  model 


Parameter 

Value 

Refs. 

f-BP 

1 x 10-3  m 

This  work 

£b 

0.8 

This  work 

£m 

0.28 

This  work 

*a,  Kc 

1.6  Wm-1  K-1 

[12,42] 

Km 

0.34  Wm_1  K-1 

[12] 

100  Wm-1  K-1 

This  work 

AL^o 

44  kJ  mol-1 

[34] 

he  au 

2672  Wm-2  K_1 

This  work 

hn2 

824  W  m-2  K-1 

This  work 

hair 

1200  Wm-2  K-1 

This  work 

Z  -  * 

li  A 

gas  gas 

Tgas* — WW — i 
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Z  =- 
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Fig.  3.  Schematic  description  of  heat  transfer  in  the  bipolar  plates  (example 
at  the  anode  side). 


flux  densities  of  the  species  come  from  the  cell  mass  transfer 
model  (Section  2.2). 

Enthalpy  transport  by  protons  and  electrons  is  neglected. 
Temperature  continuity  and  energy  conservation  at  the 
GDE/membrane  interfaces  lead  to  the  following  equations: 


97a 

Ka - 

9xa 

+  0a  — 

Xa—Pa 

dTm 

Km  Q 

dxm 

dTm 

dTc 

Km  „ 

dxm 

+  Qc 

xm=Lm 

=  —  KC - 

dxc 

with 


% — o 


xc  =0 


(47) 

(48) 


0a.c  =  e°;crpot  +  Q17  +  Q17 


(49) 


And  the  energy  conservation  at  the  GDE/bipolar  plates 
interfaces  gives: 


0a  =  -  « 


3E 


a 


a 


a*a 


xa=0 


(50) 


(51) 

0a, c  is  the  heat  flux  density  flowing  trough  the  bipolar  plate, 
which  depends  on  the  feeding  gas  and  cooling  water  temper¬ 
atures.  This  heat  flux  density  is  evaluated  thanks  to  a  resistive 
model  shown  in  Fig.  3.  The  heat  transfer  coefficients  in  the 
laminar  water  and  gas  flows  are  evaluated  with  correlation 
(52): 


0c  — 


dTc 

dXr 


(52) 


3.  Results  and  discussion 

3.1.  Overpotentials 

A  literature  review  gives  the  three  main  types  of  models 
used  to  evaluate  the  overvoltages  at  both  electrodes: 
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Fig.  4.  (a)  Influence  of  electrodes  description  on  anodic  overvoltages,  (b) 
influence  of  electrodes  description  on  cathodic  overvoltages. 

•  The  Tafel  law  results  from  experimental  observations:  this 
simple  and  widely  used  description  is  based  on  the  follow¬ 
ing  phenomenological  equation: 


Let  us  note  that  the  same  values  of  kinetic  coefficients  b[ 
and  z'o(a,c)  can  be  used  in  these  three  models  (Table  2).  Fig.  4 
depicts  the  evolutions  of  the  cathodic  and  anodic  overpoten¬ 
tials,  calculated  in  the  same  operating  conditions,  as  functions 
of  the  current  density. 

There  are  great  differences  between  the  overpotentials  pre¬ 
dicted  by  the  three  models,  which  emphasizes  the  important 
influence  of  the  active  layer  on  fuel  cell  performances.  Par¬ 
ticular  attention  should  be  paid  to  the  use  of  Tafel  law:  the 
effective  platinum  area  being  higher  than  the  ME  A  section, 
the  exchange  current  density  io  should  be  multiplied  by  a 
geometric  parameter  such  as  the  roughness  factor  y  (55).  In 
this  case  the  anodic  overvoltage  becomes  negative  (Fig.  5). 
Actually,  Tafel  law  is  adapted  only  to  high  current  density, 
far  from  thermodynamic  equilibrium 

=  bi  log  (  (55) 

\  yi 0(c,a)  J 

Some  authors  [36,37]  propose  to  change  the  value  of  Tafel 
slope  as  a  function  of  current  density,  which  has  the  advan¬ 
tage  to  be  a  simple  approach  but  does  not  provide  information 
about  the  electrochemical  processes.  On  the  other  hand, 
the  porous  and  non-porous  models  require  only  intrinsic 
parameters:  active  layer  thickness,  platinum  load  (roughness 
factor)  and  Butler- Volmer  kinetic  coefficients.  The  polariza¬ 
tion  curves  obtained  with  the  porous  and  non-porous  models 
exhibit  Tafel  behaviour.  Due  to  the  slow  diffusion  of  active 
species  in  the  homogeneous  solid  phase,  the  non-porous 
model  tends  to  overestimate  the  overpotentials.  GDE  models 
do  not  have  such  drawbacks.  Furthermore,  they  can  simulate 
the  electrode  behaviour  when  the  active  layers  porosity  is 
low  and  when  anode  hydration  is  weak  [25].  Only  this  model 
predicts  an  increasing  potential  drop  in  these  two  cases. 

3.2.  Thermal  fields 


rj  =  bi  log 


f  J cell 

V0(c,a) 


(53) 


The  numerical  results  that  are  discussed  in  this  section 
correspond  to  two  different  operating  conditions: 


•  The  models  of  non-porous  active  layer  (homogeneous 
model)  proposed  by  Jaouen  et  al.  [3]  assume  that  the  active 
layer  is  made  up  of  an  homogeneous  mixture  of  catalyst 
and  carbon  powder  embedded  in  electrolyte.  The  overpo¬ 
tential  expression  results  from  the  use  of  Butler- Volmer 
equation: 


•  In  Fig.  5a,  both  gases  are  humidified  and  preheated:  Tgas  = 
rgcas  =  80  °C,  HRa  =  0.3  and  HRC  =  0.6. 

•  In  Fig.  5b,  only  air  is  humidified  and  preheated  (Tgas  = 
70  °C,  and  HRC  =  0.9)  while  hydrogen  undergoes  no  treat¬ 
ment  (7gaas  =  25  °C,  HRa  =  0.1). 


/ 


^cell  —  P^0(a,c)  ^ 


y?  exp 


2.3)7 
bi 
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aggl  ,0  eXP 


2FD.  y 


yl 0(a,c)L  (  2.3 T) 

2FDaggly°  P  V  bi 
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/ 


(54) 


where  D^s§1  is  the  coefficient  of  diffusion  of  the  electro¬ 
active  species  in  the  homogeneous  mixture. 

•  The  GDE  models  [22,23]  take  into  account  the  existence 
of  gas  pore  (macro-pore)  and  (possibly)  of  the  ohmic  drop 
in  the  agglomerate. 


Temperature  fields  are  calculated  for  current  densities 
ranging  from  0  to  1  A  cm-2  in  the  first  case  (Fig.  5a),  and 
from  0  to  only  0.6  A  cm-2  in  the  second  case  (Fig.  5b);  this 
last  maximum  value  corresponding  to  a  cell  potential  equal 
to  zero.  The  temperature  of  the  cooling  circuit  is  set  to  80  °C 
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Fig.  5.  (a)  Temperature  fields  when  both  gases  are  preheated  and  humid¬ 
ified  (HRa  =  0.3,  HRC  =  0.6,  /Wate r  =  Tgaas  =  7gCas  =  80  °C,  /ce u  from  0  to 
1  A  cm-2),  (b)  Temperature  fields  for  dry  hydrogen  (HRa  =  0.1,  HRC  =  0.9, 
Twater  =  80  °C,  7^  =  25  °C,  T^as  =  70  °C,  /cell  from  0  to  0.6  A  cm-2). 


and  the  thermal  properties  of  the  MEA  materials  are  those 
given  in  Table  3.  The  graphs  give  temperature  distributions 
in  the  anodic  backing  layer  (0  <  x  <  230  pan),  in  the  mem¬ 
brane  (230  <  x  <  355  jam)  and  in  the  cathodic  backing  layer 
(355  <  v  <  585  pm).  In  the  following,  we  define  the  fuel  cell 
overheating  as  the  maximum  (positive)  difference  between 
the  temperature  of  the  MEA  and  that  of  the  cooling  water. 

The  linear  temperature  distributions  in  the  backing  layers 
are  characteristic  of  purely  conductive  heat  transfer.  Indeed, 
the  gas  velocities  being  low  (pfacs  a  10-3ms-1,  n{^2°  a 
10-5ms-1),  the  convective  contribution  to  heat  transport 
can  be  neglected,  as  confirmed  by  the  order  of  magnitude 
of  Peclet  number 


PH2OTm£  m  ^  t0  3 
Km 


y PjVjLm 

Ka,c 


<  1CT6 


(56) 

(57) 


Similarly,  the  temperature  distribution  in  the  membrane 
is  representative  of  conductive  heat  transfer,  but  with  a  volu¬ 
metric  heat  source. 

The  strong  dependence  of  temperature  fields  on  current 
density  is  linked  to  the  physical  mechanisms  at  the  origin  of 
the  main  heat  sources:  electrode  overpotentials  (44),  Joule 
effect  in  the  membrane  (42)  and  reaction  entropy  (43).  Nev¬ 
ertheless,  heat  generation  or  absorption  due  to  water  sorption 
or  desorption  by  the  membrane  (45)  is  proportional  to  the 
water  flux,  which  can  be  very  important  in  the  absence  of 
electric  current  when  there  is  a  strong  difference  between 
water  partial  pressures  in  hydrogen  and  air. 

The  important  difference  between  temperatures  in  the 
MEA  and  temperature  of  the  cooling  circuit  must  be  em¬ 
phasized.  When  both  gases  are  humidified  and  preheated 
(Fig.  5a),  the  temperature  in  the  membrane  and  in  the  active 
layers  can  be  a  few  degrees  above  that  of  water.  The  thermal 
model  presented  in  this  study  predicts  MEA  overheating  of 
the  same  order  of  magnitude  as  those  calculated  by  Djilali 
and  Lu  [16]  but  by  taking  into  account  the  heat  transfer  re¬ 
sistances  with  water  and  gas,  whose  temperatures  are  used  as 
parameters. 

When  the  cell  is  fed  with  dry  and  cold  hydrogen  (Fig.  5b), 
the  minimum  temperature  in  the  anodic  backing  layer  can 
be  10  °C  below  water  temperature.  This  strong  discrepancy 
is  due  to  the  cooling  effect  of  cold  hydrogen.  The  question, 
whether  such  temperature  gradients  could  induce  water  con¬ 
densation  and  thus  anode  flooding,  will  be  tackled  in  a  future 
work:  it  must  be  kept  in  mind  that  water  saturation  pressure  is 
50%  higher  at  80  °C  than  at  70  °C.  Moreover,  low  hydrogen 
supply  temperature  generates  thermal  gradients  in  the  mem¬ 
brane  (close  to  30°Cmm-1)  that  are  more  important  than 
those  due  to  Joule  effect  only  (up  to  10  °C  mm-1). 

Temperature  gradients  and  MEA  overheating  are  closely 
related  to  the  backing  layers  thermal  conductivity.  The  fibrous 
and  porous  nature  of  these  layers  does  not  allow  a  direct 
estimate  of  the  equivalent  thermal  conductivity  of  the  dry 
medium.  The  difference  between  values  coming  from  various 
sources  can  reach  one  order  of  magnitude: 

•  Gurau  et  al.  [39]  and  Wang  et  al.  [38]  used  a  equivalent 
conductivity  of  19  Wm-1  K-1  determined  with  Maxwell 
model  [40]. 

•  Kjelstrup  [41]  proposed  to  use  k^c  =  0.2  zb 
0.1Wm-1K-1,  value  estimated  from  temperature 
measurements  inside  a  single  polymer  fuel  cell. 

•  And  other  authors  [12]  agree  to  assign  a  value  of 
1.6  Wm-1  K-1  to  the  backing  layers  equivalent  thermal 
conductivity,  as  proposed  by  a  manufacturer  [42] . 

Repeated  measurements  performed  on  a  specific 
bench,  not  presented  here,  gave  us  values  of  the  order  of 
0.3Wm-1K-1  with  Quintech™  carbon  papers,  which  is 
close  to  Kjelstrup  estimation.  This  value  and  those  found  in 
the  literature  are  used  to  simulate  temperature  distributions 
in  the  MEA  in  the  same  operating  conditions  as  those  of 
Fig.  5a  (both  gas  humidified  and  preheated  and  current 
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Fig.  6.  Thermal  distributions  in  the  ME  A  for  different  values  of  ther¬ 
mal  conductivities  of  the  backing  layers  (HRa  =  0.3,  HRC  =  0.6,  rwater  = 
7gaas  =  7gcas  =  80  °C,  (a)  1  =  0.2  Wm-1  K-1;  (b)  X  =  0.3  Wm-1  K_1;  (c) 
X  =  0.7  Wm-1  K_1;  (d)  X  =  1.6  Wm-1  K_1;  (e)  X  =  19  Wm^KT1). 

density  fixed  to  0.5  A  cm-2).  The  results  are  presented  in 
Fig.  6. 

As  shown  in  Fig.  6,  the  backing  layers  thermal  conduc¬ 
tivity  has  a  great  influence  on  temperature  distributions.  The 
MEA  overheating  increases  with  decreasing  backing  layers 
thermal  conductivity,  while  the  bipolar  plates  temperature 
stays  quasi-constant  due  to  the  good  graphite  thermal 
conductivity.  A  first  conclusion  of  these  results  is  that 
the  backing  layers  thermal  conductivity  is  an  important 
parameter  to  master,  especially  because  of  the  (possible) 
presence  of  liquid  water  that  can  significantly  modify  its 
value.  Another  important  conclusion  is  that  temperature 
measurements  performed  directly  on  or  inside  the  bipolar 
plates  do  not  yield  good  estimates  of  the  MEA  temperature. 

4.  Conclusions  and  perspectives 

A  fuel  cell  model  that  takes  into  account  gas  diffusion  in 
the  porous  electrodes,  water  diffusion  and  electro-osmotic 
transport  through  the  polymeric  membrane,  and  heat  transfer 
in  the  Membrane  Electrodes  Assembly  (MEA)  and  bipolar 
plates  is  presented.  This  model  is  constructed  by  combining 
independent  descriptions  of  heat  and  mass  transfers  in  the  cell 
with  a  third  description  of  coupled  charge  and  mass  transfers 
in  the  electrodes. 

Based  on  a  Butler- Volmer  equation,  the  GDE  description 
proposed  in  this  paper  requires  only  intrinsic  parameters.  It 
allows  the  evaluation  of  the  anodic  and  cathodic  overvolt¬ 
ages  such  as  the  simulated  polarization  curves  exhibit  a  Tafel 
behaviour.  Perspectives  concerning  the  electrode  behaviour 
could  be  the  investigation  of  the  effects  of  microstructure  ge¬ 
ometry  and  of  the  active  layer  ionic  conductivity  variations, 
particularly  in  the  case  of  dry  anode.  Moreover,  modifications 
of  the  active  layer  model  will  permit  to  simulate  electrode 
flooding. 


The  results  also  show  that  thermal  gradients  in  the  MEA 
could  lead  to  thermal  stresses  at  high  current  densities.  A  dis¬ 
crepancy  between  the  air  and  hydrogen  feeding  temperature 
increases  these  thermal  stresses. 
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